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VARIATIONAL  THEORY  OF  DEFORMATIONS  OP  CURVED,  TWISTED  AND 

EXTENSIBLE  ELASTIC  RODS 

by 

Irtdj  Tadjbakhili  and  Dinitxif  C.  TjgoqdM 


ABSTRACT 

The  vaiiatlonai  theory  of  three  dimensional  deformations  of  curved  twisted  and  extensible 
elastic  rods  is  obtained  based  entirely  on  the  Idoematical  variables  of  position  and 
rotations.  The  constitutive  relations  that  define  the  reiiitive  couples  and  the  axial  force  as 
gradients  of  the  strain  energy  function  are  established.  A  candidate  for  the  strain  energy 
function  on  the  basis  of  classical  assumptions  is  presented. 

INTRODUCTION 


In  the  Historical  Introduction  of  the  "A  Treatise  on  the  Mathematical  Theory  of 
Elasticity, Love  (1893)  narrates  that  in  1742  Daniel  Bemonlli  wrote  to  Euler  suggesting 
that  the  differential  equation  of  the  elastica  could  be  found  by  making  the  integral  of  the 
work  done  or  the  square  of  the  curvature  a  minifflum.  Acting  on  this  suggestion  Euler  was 
able  to  obtain  the  differential  equation  of  the  elastica  and  the  various  forms  of  it.  Thus  the 
concept  of  the  strain  energy  was  bom  and  the  foundations  of  the  variational  theory  of 


Love  are  applicable  to  an  initially  bent  and  twisted  rod. 
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Our  aim  in  this  paper  is  to  establish  a  variational  formulation  for  the  title  problem  and  in 
the  process  infer  the  existence  of  the  strain  energy  function  and  determine  the  constitutive 
relations  that  relate  this  function  to  the  bending  and  twisting  couples  and  the  axial  force 
within  the  rod.  This  development  together  with  equations  of  equilibrium  and  the  geometry 
of  deformation  define  a  direct  approach  and  an  exact  nonlinear  theory  for  the  three 
dimensional  deformation  of  a  one  dimensional  elastic  medium  capable  of  resisting  bending 
twisting  and  extension.  Going  a  step  further,  in  order  to  actually  construct  an  explicit 
form  for  the  strain  energy  function,  we  enter  the  realm  of  hypothesis  and  use  KirchhofTs 
description  of  deformations  in  a  thin  rod.  Thii  view  enables  ui  to  determine  a  strain 
energy  function  that  can  be  used  in  engineering  applications. 

The  recent  history  of  investigations  of  the  rod  theories  consists  of  developments  along  two 
separate  streams,  the  direct  approach  and  approximations  from  three-dimensional 
continuum.  In  the  direct  approach  a  one-dimensional  continuum  view  is  pursued  and  the 
medium  is  supposed  endowed,  in  addition  to  its  position,  with  vector  fields,  the  directors, 
that  are  to  be  interpreted  appropriately  to  define  bending,  twist  and  extension  properties  of 
a  rod.  This  approach  has  its  origins  in  the  work  of  E.  and  F.  Cosserat  (1909)  and 
numerous  investigators  have  contributed  to  it  among  them  Naghdi  (1983),  Naghdi  and 
Rubin  (1984),  Whitman  and  DeSilva  (1970),  Green  and  Laws  (1966),  and  Eriksen  (1970). 
Extensive  investigation  into  the  qualitative  aspects  of  the  nonlinear  theory  such  as 
questions  of  existence  of  solutions  and  global  behavior  have  been  carried  out  by  Antman 
(1976).  His  basic  work  entitled  "The  Theory  of  Rods"  (1973)  describes  these  theories  both 
as  approximations  to  the  three-dimensional  continuum  theory  and  as  a  one-dimensional 
continuum  with  directors.  The  work  presented  here,  although  pertains  to  a 
one-dimensional  continuum,  does  not  use  directors,  but  is  formulated  entirely  on  the  basis 
of  Idnematical  quantities  eonsilting  of  the  positioa  vector  of  points  along  the  curve  of 
centroids  and  the  orientation  angles  of  the  cross  sections  of  the  rod  relative  to  a  fixed 
coordinate  system.  It  is  a  generalisation  of  the  work  of  Tadjbakhsh  (1966)  in  which  the 
theory  of  planar  motion  of  the  extensible  elastiea  wu  described. 

The  history  of  construction  of  approximate  theories  in  the  context  of  three-dimensional 
nonlinear  continuum  theory  is  also  varied  and  to  it  many  investigators  including  some  of 
the  above  authors  have  contributed,  see  for  example  Green,  Naghdi  and  Wenner  (1974). 
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KINEMATICS 

An  elastica  is  a  nearly  uniform,  slender  rod  of  finite  length.  In  the  unstressed  state  the 
centroids  of  the  cross  section  form  a  space  curve  C  that  is  called  the  reference  curve  with 
an  arc  length  s.  The  orientation  of  the  principal  axes  of  the  cross  sections  vary 
continuously  along  the  rod.  This  means  that  in  the  unstressed  state  the  rod  has  arbitrary 
twist  and  curvatures.  With  respect  to  a  fixed  Cartesian  frame  the  coordinates  of  points  on 
the  reference  curve  are  denoted  by  Xi(s),  i  «  1,  2,  3.  The  cross  sectional  area  can  be  a 
slowly  varying  function  of  s  and  will  be  denoted  by  A(s). 

Attached  to  any  point  s  of  C  a  Cartesian  coordinate  system  with  axes  yi,  i«  1,2,3  will  be 
assumed  and  will  be  referred  to  u  the  body  axes.  As  shown  in  Fig.  1  the  yi  axis  will  be 
pointing  in  the  direction  of  the  increasing  s  and  yi  and  ys  axes  will  point  along  the 
principal  axes  of  the  cross  section  such  that  yi  -  frame  is  a  right  handed  system. 

As  the  rod  deforms  the  curve  C  acquires  new  configuration  c.  The  arc  length  along  c  is 
denoted  by  ^(s)  to  account  for  extensibility  of  the  rod.  The  position  of  a  point  s  of  C  that 
in  the  unstressed  state  is  X(s)  at  the  stressed  state  will  be  x(s),  both  referred  to  the  same 
fixed  Cartesian  frame.  Noting  that 

('4)  (1) 

tlu  ittain  e  ii  dcflud  bf 


*  =  {'-1  (J) 

t 

where  e>0  denotes  extension  and  e<0  contraction.  The  strict  positivity  of  (  implies 
-1  <e<». 


Two  sets  of  orthogonal  unit  vectors  e i  and  ni  will  be  assumed  for  the  body  frame  y  and  the 
fixed  frame  x  respectively.  Denoting  by  l||(s)  the  elements  of  the  matrix  of  direction 
cosines  between  these  two  sets  of  unit  vectors  one  has 
ai  •  Iqij 


(3) 
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Ci  »  Ijl  &J  (4) 

Let  us  consider  an  element  of  the  rod  Ads,  Fig.  I,  and  use  the  dynamical  analogy  of  £.1. 
Routh.  Love  has  noted  that  if  the  frame  were  to  move  with  unit  speed  along  the  curve  c 
such  that  at  any  point  ^  of  c  it  has  the  orientation  of  the  j  frame  at  that  point  then  the 
angular  velocities  ui  and  uit,  about  the  yi  and  yt  axes  respectively,  will  be  the  principal 
curvatures  /ti  and  /<]  of  the  rod.  Also  the  angular  velocity  a;s  will  be  the  twist  curvature  Ki 
of  the  tod.  Thus  the  formulas  that  define  angular  vdodties  from  direction  cosines  can  be 
used  to  determine  curvatures,  provided  time  differentiation  is  replaced  by  differentiation 
with  respect  to  t  The  curvature  vector  is  defizmd  by 

«  ■  ei  (5) 

with 

«i=>7i8h^ljh  (®) 

where  =  etjk  (eijk  +  l)/2,  (no  sum  on  i,  j,  k),  and  stjk  is  the  alternator  tensor  with  its 
non-eero  components  being  tias  •  ■  «iai  ■  +1  wd  «iia  ■  <iai  ■  ■  “I-  Noting 

that  d/d(  »  (1  +  d/ds,  we  define  the  three  refisrence  curvature  parameters  ki  by 

ki  *  (l+e)Kt  (7) 

For  future  use  we  record  the  formulas  Ibr  derivatives  of  the  unit  vectors  ei 

e\  m  eiiji  kj  ek  (8) 

Associated  with  the  unstrained  state  of  the  rod  there  will  be  the  unit  vectors  Bi  that  bear 
the  same  relationship  to  the  principal  directions  of  the  cross  sections  of  the  rod  as  do  ei  in 
the  strained  state.  In  particular,  if  Ly  are  the  direction  cosines  of  Ei  with  respect  to  m  and 
Ki  are  the  curvatures  and  twist  of  the  rod  ia  the  unstrained  state,  then  formulas  (3)  (7) 
apply  with  lij,  et  and  k{  replaced  by  L^,  Ei  and  Ki,  reipectivtiy,  and  e»0. 

Since  es  and  are  both  unit  tangents  to  the  ^tral  line  one  has 
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Xj  a  (l+e)lij 

whfre  xj  represent  the  components  of  x  in  the  fixed  fr&me  i.e. 

XaXiHi  (10) 

Thus  (9)  is  a  differential  constraint  that  relates  the  coordinates  of  the  central  line  to 
direction  cosines  of  the  principal  axes  of  the  cross  section  of  the  rod. 

The  strain  e  and  the  three  curvature  parameters  ki  form  a  set  of  four  field  variables  that 
completely  characterises  the  deformation  of  the  rod.  In  subsequent  sections  conjugate 
forces  will  be  defined  through  the  standard  procedure  of  differentiating  the  strain  energy 
fonction  with  respect  to  these  kinematieal  vaiiaUes. 


EQUATIONS  OF  EQUIlIBBIUM 

Referring  to  the  body  set  of  axes  ei  one  can  define  the  vector  F  of  the  resultant  shear 
stresses  F  i  and  Fj  of  the  cross  section  and  the  axial  stress  resultant  Fs.  Similarly  one  may 
define  the  resultant  couple  stress  vector  If  consisting  of  the  bending  moments  Mi  and  Ma 
about  the  axes  yt  and  ya  and  the  torque  Mf.  Explicitly  we  have 


F-F,ei 


(11) 


and 


y»Miei  (13) 

The  well  known  force  and  moment  equilibrium  equations  are 

P^+f-0  (18) 


and 


(W) 
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where  f  and  m  repreaent  the  distributed  force  and  couple  acting  on  the  rod,  per  unit 
undeformed  length  s. 

The  scalar  components  of  equations  of  equilibrium  can  be  expressed  in  the  body  set  of  axes. 
For  this  purpose  one  needs  to  express  all  vector  quantities  in  terms  of  unit  vectors  ei  and 
use  (8).  Then  (13)  becomes 


F  j  +  k jF  s  —  k|F 3  +  f^  a  0 

(Ua) 

F,  +  k|Fi-kiF|  -I-  fjs  0 

(IJb) 

F|  +  kiP3  —  ksF j  +  fj  »  0 

(I5c) 

while  (14)  assumes  the  form 

M  j  +  kaM|  —  k|Ma  ■*  (1  +  e)F  3  +  m^  ■  0 

(16a) 

m'  +  ksMi  -kiMi  +  (1  4-  e)Fi  4-  m][  -  0 

(16b) 

M|  ^  k^Ma  —  kaMi  ^  m|[  *  0 

(16c) 

Here  the  superscript  7  on  the  components  of  f  and  m  denote  these  components  in  the  body 
reference  frame. 

To  express  the  equations  of  equilibrium  in  the  fixed  frame  we  introduce  the  components  of 
the  stress  resultants  in  that  frame.  Thus 

fJ-IijP,  ,  x’-IgM] 


(IT) 
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Then  (13)  becomes 

Fr+ft  =  0  (18) 

and  (14)  assumes  the  fonn 

Mf  -  (l+e)(lia  F*)  +  m*  ■  0  (19a) 

M*'  +  (l+e)(lii  Fi)  +  m*  »  0  (19b) 

MJ'  +  mj  »  0  (19c) 


CONSTITUTIVE  RELATIONS 

Either  one  of  the  set  of  equations  (15)-(16)  or  (18)-{19)  can  be  considered  as  the  governing 
differential  equations  of  equilibrium  of  the  rod>  These  equations  have  to  be  supplemented 
with  constitutive  relations  that  define  resultant  axial  stress  Fi  and  the  resultant  bending 
and  twisting  couples  Mi,  Mj,  Ms  in  terms  of  the  axial  strain  e  and  curvatures  ki,  k2,  k|. 
For  this  purpose  we  assume  the  existence  of  a  strain  energ7  function  W(e,  kt)  per  unit 
undeformed  length  which  is  invariant  under  rigid  body  translations  and  rotations  in  the 
deformed  configuration.  We  consider  the  variational  problem  of  minimum  of  the  total 
potential  energj  u  the  equivalent  of  equilibrium  states  of  the  rod.  This  potential  can  be 
expressed  as 


*(.  •)  -  /  V(*.  ki)  +  [*i  - 

0 

,  _  -  -  -  -  L 

—  fj  Xj  —  mj  iP(}ds  [Mj  ^  Fj  xj  ^  (20) 

where  f|  ,  mJ  and  are  the  applied  forces  and  moments  at  the  end  s  ■  0  and  s  «  L 

respectivelF. 
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The  functions  Xi  are  Lagrange  multipliers  that  allow  the  constraint  (9)  to  be  incorporated 
within  the  functional  U.  As  a  result  xi  and  lij  can  be  regarded  as  independent  variables. 
Additionally,  the  constraint  (9)  implies  the  definition  (1)‘^2)  for  the  strain  e  and  hence  in 
(20)  e  can  also  be  viewed  as  an  independent  variable.  To  see  this  we  need  to  note  that  if 
each  side  of  (9)  is  multiplied  by  itself  we  obtain  =  (l+e)2lijlij  =  (1+e)^  which  is  a 

restatement  of  (l)-{2).  The  terms  f^i  and  in  the  integrand  of  (20)  represent  the 
density  of  the  potential  of  the  applied  forc^  and  moments  on  the  rod.  The  angles 
represent  the  rotations  from  Ei  to  ei  when  these  directions  are  assumed  to  issue  from  a 
common  origin.  These  angles  are  determined  through 

cos^i  *  ei  •  Ei » Iji  Lji ,  i»l,2,3  (sum  only  on  j)  (21) 

The  direction  coslns  Ijt  are  characterised  by  three  orientation  angles  Bt  that  can  be 
selected  in  a  variety  of  ways  and  represent  three  finite  rotations  about  the  unit  vectors  ei 
or  m.  If  these  rotations  are  properly  selected  any  initial  orientation  of  a  cross  section  may 
be  brought  to  any  arbitrary  final  orientation.  Kane  et  al.  (1983)  list  at  least  24 
possibilities  for  order  of  rotations  of  the  anglea  $i,  h  Bt  about  the  body  set  of  unit  vectors 
ei  or  the  fixed  set  of  unit  vectors  m.  Thus  in  carrying  out  variation  with  respect  to  a 
particular  <fii  (say  ^t)  we  need  to  select  a  particular  sequence  of  rotation  angles  the  last 
of  which  coincides  with  Same  procedure  must  be  used  when  carrying  out  variation  with 
respect  to  or  tp^. 

With  these  preliminaries  we  note  that  the  Euler  equation  corresponding  to  variations  dxi  is 

simply  A(  -h  f^  ■  0,  which  when  compared  with  (16)  ideatifiei  At  with  F*  Next 
considering  the  variation  with  respect  to  e  we  obtain 

(22) 

which  is  the  constitutive  relationship  determining  the  axial  force  Ft  as  the  derivative  of 
strain  energy  with  respect  to  axial  strain  e. 
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We  now  turn  to  the  Euler  equation  corresponding  to  the  variation  Stfii.  For  orientation 
angles  of  the  cross  section  we  select  the  sequence  of  body  rotations  first  ffjej,  second  ^^es 
and  third  »  (fiiCi-  Then  the  matrix  1  of  direction  cosines  is  given  by 


CjC|  -^iCjSi+SiSa  SiCjSs+CiS] 
S3  C1C3  “SiCi 

SaC|  CtSaSi+SiCs  — SiSjSs+CiCj, 


where 


(23) 


Ciicoidi  ,  Siisln^i 


'10  0 

'  cold  0  lind 

Mfl‘ 

0  cos^  -Sind 

.  B(*)« 

0  1  0 

0  Sind  cold 

^-llnd  0  00s d 

‘cosd  -Sind  O' 

C(tf)- 

Sind  cold  0 

0  0  1 

*  • 

Subsequently  we  find  from  (7) 
S|  + 

k, «  $2  +  ^|8i 

k|  s  -  SjC|  +  ^iCi 

From  (20)  we  have 


(24) 


(25) 


(26a) 

(26b) 

(26c) 


0 


(2T) 
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Noting  that  dfd^i  s  d/d8i,  we  have  from  (28)  dk\Jd8\  =  0,  =  ks,  dk.il *  '■1^3 

and  from  (23)  dlis/d^t »  Using  alio  the  inverse  of  (17a),  (27)  becomei 


+  -  (1+«)F)  + 


(28) 


In  exactly  the  same  manner  one  can  proceed  to  determine  the  Eulex  equation  for  a 
variation  Now  the  consecutive  sequence  of  body  rotations  8303,  9\ti,  is  selected 
with  9i  m  ipi.  Without  going  into  details  we  obtain 


(29) 


For  variation  of  we  adopt  the  sequence  of  body  rotations  ^lei,  ^aea,  ^ler  with  $t  »  <fit- 
The  matrix  of  direction  cosines  is 


l-A(^i)B(^3)C(  Go¬ 


eses 

SiS  jC  j+CiSi 
.-CjSaC  j+SiSi 


-C3S1 

-SiSaSi+C|C| 

C|SaSs+SiC| 


Sj  * 
•S 1C3 
CIC3, 


(30) 


with  curvatures  given  by 

k(  ■  liCaCa  +  l|S| 

(31.) 

k)  ■  #1  Ci  -  l|CaS| 

(31b) 

k|  *  l|  +  IjSa 

(31C) 

For  this  case  lu  does  not  depend  on  la  and  hence  the  Enler  variational  equation  assumes 
the  form 


(32) 
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Comparison  of  equations  (23),  (29)  and  (32) 
establishes  the  constitutive  relations 


with  equations  (I6a)  -  (16c)  respectively 


The  specified  boundary  conditions  must  be  consistent  with 

((f  i  -  F!)&Ci  +  (101  ^  .  0 

which  imply  that  at  a  boundary  point  either 

Fi-Fi  or  xi«xi 
and 


(33) 


(34) 


(33) 


Mi  a  101  or 


where  x,  and  axe  specified  position  and  orientation. 


(36) 


A  STRAIN  ENERGY  FUNCTION 

T  “•"*  »“»  *“«•«  -  “"-iT  tk. 

Ki«^  t ^  O'  «“  "Xl-  for  tto  ,„P0H  «e  i.roke  th, 

*  V  «!»»«»  pUm  aoM  Hctiou  of  rod  .U»  u.  norni.1  to 

Mfort  tad  tile  deforoaUoa  eta  be  |fm  by 


a-X(t)  +  y,E,(i)  +  y*(,) 
tad 


(37) 


f-i(«)  +  <»(l)lyi»,(,) +  ,,,,(,)) 


(33) 
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respectively.  The  parameter  a((8)  is  to  be  fixed  by  enfordng  tractioa-free  boundary 
conditions  on  the  lateral  surface  of  the  rod. 

I’sing  the  concept  of  extensional  strains  for  stretching  of  line  elements  and  distortion  of 
angles  between  perpendicular  lines  as  shear  strains  (Wempner,  1991),  we  define 
components  of  strain  by 

*  i(si  *  fii  “  Gi  •  Gj) 

where 

^  Sa  -  ^  •  «a. 

dr 

^  *  ^'7 fix  +  a'ffii  +  +  +  (l  +  e)as  (40) 

^  Ga  ■  ^  •  Ea.  G|  «  ^  -  j{r[  +  yjBj  +  E|  (41) 
Using  (8)  we  can  establish 


®i  '  ®j  ■ 

(42) 

•  *1  ■  <ij A 

(43) 

where  Ki  is  the  eurvatuei  and  twist  in  the  unstrained  state.  Therefore  (9)  yields  as  strain 
components 

*u  ■  <11  ■  ♦  («*  ”  1) »  «ij  ■  0 

«ii  -  ♦[«(«'7|-ay|k,)  +  y,KJ 
<11  ■  (o'yi  +  «  yaks)  -  yjK  J 

<is  -  e  +  ie*  -  ytl(l+e)a  ka  -  Ka)  +  ya((l+e)o  k,  -  K  J  -  yja  (nnnki  -  KiKa) 

+  iySa'  *  +  a*  (k}  +  ki)  -  K  j  -  K})  +  ^y*  [o'  a  +  oa(k}  +  kj)  -  K{  -  Kj  (44) 


01/14/93  15  25  0409  845  6051 


TAMU  Aerosp.  Eng. 


il014/018 


13 


For  a  linear  isotropic  elastic  material  the  non-eero  stress  components  per  unit  strained 
area  are 


ffU  s  (A+G)  (aM)  +  Xt»  , 

<y»a-i  0 

ffii  s  (A+G)(<»a-1)  +  A««  , 

»  2G«,t 

-  A(aM)  +  (A+2G)  <« , 

»  2GS|i 

(45) 

where  X  and  G  are  Lame's  conatant  and  the  shear  modnlns,  respectively.  The  traction  per 
unit  undefonned  area  is  then  ^ven  by 

-  j«fi  (46) 

One  can  define  the  asdal  stress  resultant  F]  by 

F|  -  dA  - 

A 

*  A(l-fe)  ^(A+2G)(e+^*)  +  A(o^l)j 
+  {X-t-id)  [t.((l.M)akrKJ(iki  + 1,1(1  4-  •)(ik,-KJak, 

+  ^  I,  (1  4- t)(a' >  4- <A1  -  K*  4- (A;  -  KD 

+  iI,(1  +  «)(«''4.<A|-K|4.  A}-K{)]  (47) 

where  Ii  /^y|  dA  and  It  •  /^y’  dA.  Similarly  we  have 

M,-  dA  >  a(A4-90)  lj(l  4>  «)(*,-  K,j(l  4- 1) 

4-  1,iA,[a(«>-1)  4-  (A4.J<J)(44-4«>)]  (4*») 

M,  -  -^a  7,t'.*,dA  -  a(A  4-  JO)  I,[(l  4-  #)i*,  -  K,]  (1 4- «) 

4-  I,a>k,[A(a»-l)  4-  (A4.2Q)(e4.H’)] 


(41b) 
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M,  -  dA  -  JGa^CcAj-Kj) 

+  4-  (A+2G)(e+4e2)]  (48c) 

where  One  may  note  that  the  integrability  conditions 


are  satisfied.  Hence  existence  of  a  strain  energy  function  is  assured  and  by  integration  we 
have 

W  «  A  e*(l+e+|^)  -  AA  (l-aJ)(e  +  ie*) 

+  A  I,a»  (a*-l)  ^  +  (A+3G)  Ii  a  k,  [l±2|2±lS  ok^  -  (1+«)K,] 

+  A  IjO*  (cf^l)  ^  +  (A+2G)  I|  a  kj  ~  (l+e)Kjj 

+  (1  +  e)*[li  (a'*  +  cA?  -  K5  +  Ij  (o'«  +  (A|  -  K|)] 

A  J  (oM)  J  (e+*ea)  ((A}  -  K|)  +  J  G  Ja  (ok,  -  K,)* 

(50) 

For  an  initially  straight  rod  Ki  should  be  set  equal  to  aero  in  (47),  (48)  and  (50).  The 
above  form  of  W  reflects  matetial  isotropy,  i.e.,  W(e,kiJ[i,ks)aW(e,k|,kuki),  provided  that 
Ii-Is. 


We  note  that  poaitiTe  curraturei  imply  positive  bending  moments  and  oonversevtfy 
negative  cumtuies  imply  negative  bending  moments  provided  that  e>(-l+l/^)  for 
(k  •  1.  This  shows  that  equations  (48)  have  a  limited  range  of  validity  if  the  sense 
correspondence  between  moments  and  curvatures  is  be  to  retained.  We  also  note  the 
second  order  coupling  between  the  squares  of  the  eurvatuns  and  the  axial  strain  e  in  (47), 
which  implies  that  axial  force  can  be  generated  by  bending  or  twist  only. 
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The  parameter  a(s)  depends  on  the  boundary  conditions  applied  at  the  lateral  surface  of 
the  cod.  If  the  lateral  surface  is  fixed,  then  a(s)  -  I.  For  zero  tractions  on  the  lateral 
surface,  a  condition  appropriate  for  thin  flexible  rods  is  adopted  according  to  which  the 
average  of  and  cr>3  over  the  cross-section  should  vanish,  i.e. 

a  j[ff*SdA  s  0  (51) 

The  above  condition  reduces  to 

H(a,  k,.  e)  s  a' J  +  a»(I,k!  +  I,ki  +  ft}  +  A) 

-  (I,K}  +  I,K1  +  JKj)  +  2Ae{l  +  ♦«)  -  A  -  0  (52) 

This  equation  should  be  interpreted  as  a  di&rential  equation  for  a(s)  when  ki  and  e  are 
known.  To  achieve  this  (52)  is  solved  in  an  iterative  procedure  in  which  at  every  step  k} 

and  e  are  known.  We  begin  by  writing  (52)  as  li^^  ,  kf  ,  eg)  «  0  with  n  *  0, 1,  2, 

3, ....  For  the  first  iteration  (a  ■  0),  e, «  0  ,  kj «  K^,  otj  ■  1,  and  the  six  equations  (15)  - 

(16)  after  using  (22),  (33)  and  (50)  contain  only  six  unknown  quantities  F|,  Fj,  ei  and  k| 
when  the  externally  applied  forces  and  moments  are  prescribed.  Solution  of  this  set  of 

equations  enables  one  to  use  k{  in  the  curvature-orientation  angle  relations  such  u  (31)  to 
determine  the  latter  i.e.  Now  (^)  are  known  and  one  proceeds  to  determine  from 
(21)  and  X|  from  (9)  using  the  appropiitte  boundary  conditioas.  The  solntion  for  the  first 

iteration  is  complete.  One  eaters  (52)  with  ft  and  k|  and  oomputei  o,  and  the  iteration 
proceeds. 
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